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Abstract 

Electrically generated spin accumulation due to the spin Hall effect is imaged in n-GaAs channels 
using Kerr rotation microscopy, focusing on its spatial distribution and time-averaged behavior in 
a magnetic field. Spatially-resolved imaging reveals that spin accumulation observed in transverse 
arms develops due to longitudinal drift of spin polarization produced at the sample boundaries. 
One- and two-dimensional drift-diffusion modeling is used to explain these features, providing a 
more complete understanding of observations of spin accumulation and the spin Hall effect. 
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The observation of current induced spin polarizatioiP and the spin Hall effectP demon- 
strate all-electrical generation of spin polarization through the spin-orbit interaction. A 
detailed understanding of spin accumulation produced by these phenomena is a prerequisite 
to advancing spin-based electronics applications. Sample boundaries in spin-orbit systems 
are an active topic of research, affecting both the observed spin accumulation near and 
away from sample edges 3 and even the proper theoretical definition of a spin current 4 . For 
spin Hall systems in the extrinsic regime, many of these complications can potentially be 
avoided because the spin-orbit coupling is weak, enabling simple solutions to spin diffusion 
equationP which have been used to explain spin accumulation in a number of studie^^. 

Despite these expected simplifications in the analysis of the extrinsic spin Hall effect, 
spatially-resolved optical measurements have revealed spin accumulation which is not en- 
tirely explained by simple one-dimensional spin diffusion models 2 3 . In this Letter we measure 
spin accumulation with unambiguous contributions from both drift along electric fields and 
spin precession about an applied magnetic field. We develop a diffusion model in one- and 
two-dimensions which can reproduce the observations as well as allow reinterpretation of the 
details of recent experiments on the spin Hall effect. These results provide a more complete 
understanding of the behavior of electrically generated spin polarization in micron-scale 
devices. 

The structures discussed in this manuscript are fabricated from a 2-/xm thick epilayer of 
n-GaAs deposited on 200 nm of undoped Alo.4Gao.6As grown by molecular beam epitaxy. 
The n-GaAs is Si-doped with n = 3 x 10 16 cm -3 (at 30 K) for long spin coherence times 
as in Ref. El Using standard photolithography and wet etching, the n-GaAs is patterned 
into 300-/xm long, 60-/xm wide channels with a 60-/im long, 30- /im wide transverse arm (Fig. 
la). The n-GaAs is contacted by annealed Au/Ni/Au/Ge/Ni. The devices are mounted in 
a helium flow cryostat with the channel defined to be the x direction and a magnetic field B 
applied in the y direction. All measurements reported are at 30 K, and have been repeated 
in a variety of similar n-GaAs wafers and channels of various dimensions. An AC square 
voltage with frequency fy = 1.168 kHz is applied across the device, driving a current in the 
±x direction and generating spin polarization s z in the ±z direction due to the spin Hall 
effect 2 . 

Spin accumulation is measured using a Kerr rotation microscope with a spatial resolution 
of rsj 1 /am 8 . A linearly polarized beam from a pulsed Ti: sapphire laser tuned to A = 821 nm 
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is focused through a microscope objective normal to the sample surface and the polarization 
of the reflected beam is analyzed with a balanced photodiode bridge and lock-in amplifier 
operating at fy. The axis of linear polarization rotates by an angle 6k proportional to s z 
due to the Kerr effect. The microscope objective is scanned over the device to spatially map 
s z while using the reflectivity of the laser over the channel edges to calibrate the coordinate 
position. 

Spatial profiles across the channel width (Fig. lb) show a hyperbolic sine form for s z (y) ) 
as expected from solution of basic diffusion equation^El. Spin accumulation in the channel 
is measured at least a few spin diffusion lengths (L s ~ 5/xm) away from the arm edge to 
minimize diffusive effects from the transverse arm. Sweeps of B produce a quasi-Lorentzian 
lineshape due to the time-averaged projection of s z known as the Hanle effect 2 (Fig. lc). 
The spin coherence time can be extracted from the half- width Bi/ 2 of a Lorentzian fit to the 
magnetic field scans, r = hj ' g\iB^\i2^ where h is the reduced Planck's constant, g = —0.44 
is the g factor for GaAs, and \ib is the Bohr magneton, r decreases with increasing applied 

which has been previously explained as increased D'yakonov-Perel spin scattering due 
to electron heating 9 . As a function of y, the central peak narrows away from the channel 
edge, causing an apparent increase in the extracted r toward the center of the channel 
(Fig. lb inset). This feature is consistent with previous observations of the spin Hall 
effect 27 . Increased averaging in the current measurements also reveals that s z changes sign 
away from the central peak. Similarly, magnetic field scans in the transverse arm are more 
complicated than those in the main channel, exhibiting sign reversals and narrow central 
peaks (Fig. Id) that cannot be fit by a Lorentzian lineshape. The cause of the departures 
from simple Lorentzian behavior in both the channel and the arms is expected to be due to 
a combination of spin precession, drift, and diffusion^, but detailed calculations of both the 
spatial and B dependences of the experimental data have so far been lacking. The remainder 
of this Letter will present modeling necessary to understand the above observations of spin 
accumulation due to the spin Hall effect. 

In order to model spin accumulation, we assume that spin is generally conserved except 
for a decay with time constant r. The sample boundaries are treated as hard walls; neither 
charge nor spin can penetrate sample edges, and the edge causes no spin-orbit induced spin- 
flips. These approximations are both reasonable as GaAs is expected to be in the extrinsic 
limit of the SHE where spin-orbit coupling is weak. Under these assumptions, the vector 
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spin density s l obeys a steady state continuity equation with additional spin precession and 
decay terms: 

_ = -v.J!- 7 + ( i JxU- = o (i) 

where uj = gjiBB/h is the Larmor frequency vector and J l s is the spin current of the z-th 
component of spin, which is defined in the simplest fashion as a sum over all particles of the 
product of spin and velocity: J\ = s l v. A general treatment using a more complicated 
definition of spin currentF, which is potentially necessary in the intrinsic regime, is not 
necessary in this case. For electrons with isotropic, spin-independent transport properties, 
the spin current tensor is: 

Jl d = -(13% - Ddjj + a SH e^ k E k (2) 

where \± is the charge mobility, E the electric field, D the diffusion constant, gsh the spin 
Hall conductivity, and e^ k the antisymmetric symbol. The first term represents drift of spin 
polarized carriers along the electric field, the second term governs diffusion of accumulated 
spin polarization, and the third is the spin current generated transverse to an electric field 
by the spin Hall effect. 

Spin and charge coupling 5 is assumed to be negligible for simplicity, decoupling the electric 
field solution from the spin accumulation. This additional coupling does not appear to be 
necessary to capture the salient features of the experimental observations. Introducing the 
spin diffusion length L s = \J Dt^ the drift-diffusion equation in the channel is obtained from 
Eq. [l]and Eq. [2} 

aV + ^w-£ + ^x^=o (3) 

In Eq. 3, we assume the system is in the linear transport regime, where osh is independent 
of electric field and the electron density is uniform. Because the electric field has no curl, the 
divergence of the spin Hall current is identically zero: dj(asH^ k Ek) = (JsH{^ k djEk) = 0. 
For linear transport, the spin Hall current does not act as a source for spin accumulation 
in bulk; all spin accumulation arises from boundary accumulation at the channel edge. Re- 
cent modeling' has shown that bulk spin accumulation can indeed occur in a non-linear 
transport regime where the divergence of the spin Hall current is not identically zero; we do 
not study this effect in this manuscript. Eq. [3] represents a system of coupled linear differ- 
ential equations for s l written entirely in terms of experimentally determined parameters. 
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These equations include drift and spin precession terms, both of which have demonstrable 
consequences in the experimental data presented in Fig. 1. 

In the case of an infinitely long channel oriented along the x direction (which is applica- 
ble far from the contacts and transverse arms), s l will be independent of x and the local 
electric field will be a constant E = Ex, causing the drift term to vanish and validating the 
assumption gsh{E) = gsh of Eq. [3j Specializing to the case of B along y, Eq. [3] reduces 
to two non-trivial coupled differential equations for s x (y) and s z {y) ) which can be solved 
exactly by enforcing the boundary condition h • J\ — 0. These diffusion equations including 
spin precession were solved for the case of spin accumulation in a one-dimensional system 
with polarized spins injected from ferromagnetic contacts^. The difference in this case is 
that the spin polarization is generated at the channel edge by the spin Hall effect. The 
spatial profile s z (y) at B = is the same hyperbolic sine solution used in earlier analysis 
of the SHE 5? . Near the channel edge, s z (B) has the expected Lorentzian shape from the 
Hanle model. Away from the channel edge, spin precession modifies the steady state solution 
s z (B) ) resulting in a narrowing central peak and oscillations in the tails. Calculations of the 
s z (B) dependence are shown in Fig. 2a, assuming r = 2.5 ns and L s = 4.5/xm, overlaid with 
the normalized data from Fig. lb. The narrowing of the central peak and the sign change 
of s z are evident both in these B scans and and the full B and y dependence of the data 
and calculations in Fig. 2c. Within the diffusion model, the central peak narrowing and 
the non-Lorentzian tails of Fig. lc can be entirely understood as spin precession of diffusing 
spins with a constant r without invoking a position-dependent spin coherence time. 

We apply the model discussed above to two-dimensional geometries as well, focusing on 
the transverse arm originally investigated by Sih et al. 3 In this geometry, the drift term of 
the equations plays an important role. Two-dimensional finite element modeling is used to 
first solve for the electric field E(x,y) and then use Eqs. 1 and 2 to find the magnitude of 
spin polarization s z (x ) y). As before, we fix L s = \JDt to be equal to that measured from 
magnetic field scans along y in Fig. lb. Since the magnitude of the electric field, E : varies 
throughout the sample, we use in the model the r(E) and /jl(E) measured separately for this 
sample. We note that good agreement with the data is also observed by fixing r and \i to 
their E — values, with the reduced mobility being countered by the factor of five increase 
in r. The results we present here use the measured r(E) and /jl(E), however, as this more 
likely describes the true dynamics of the mobile electron spins. The model accounts for both 
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the B = spatial profile and the field dependence observed in linecuts at x = (Fig. 2d). 
Contrary to the model in Ref. 3 ; this model does not require electron spins to have a drift 
velocity transverse to the electron current. 

In order to confirm the effects of drift in accumulation from the spin Hall effect, we 
obtain two-dimensional maps of spin polarization in the transverse arm by modulating the 
magnetic field at = 3.3 Hz and detecting the resulting Kerr rotation with lockin amplifiers 
at fv ± 2/b. The resulting signal will be proportional to the second derivative of the 
magnetic field profile near B = ? Kato:2004b), allowing efficient two-dimensional imaging 
of spin accumulation without time-consuming magnetic field scans. Fig. 3a shows an image 
of KR in the 60- /im long transverse arm extending from a 300-/xm long channel. Spin 
accumulation is observed across the entire width of the arm, extending 40 /im from the 
main channel, as observed in Ref. [31 This measurement is performed with a symmetric 
square wave current, so the observed accumulation is symmetric about the center of the arm 
at x = (Fig. 3a). Conclusive evidence of the importance of drift in the spin accumulation 
is obtained by applying an AC voltage with a DC offset to create a square wave electric field 
between mV//xm and E. This slight modification allows the lock-in detection technique to 
remain unchanged, but removes the symmetric current flow so that electrons only flow in one 
general direction along the channel. Figs. 3b and 3c show the data and model for the spin 
accumulation in the transverse arm under unidirectional current flow. The accumulation 
profile in the transverse arm is asymmetric about the arm center, and reverses direction 
when the current direction is reversed. The spin polarization observed in the transverse arm 
is neither generated in the arm itself, nor does it come from the spin Hall current generated 
in the bulk of the channel; it consists primarily of spin polarization sourced near the channel 
edge that travels longitudinally along the electric field lines that penetrate into the arm. 
These results accentuate the importance of including spin drift along the electric field in 
analysis of both one- and two-dimensional spin accumulation profiles. 

In conclusion, we have provided experimental evidence that, in addition to diffusion, 
drift and spin precession are important factors in understanding the spin polarization gen- 
erated by the spin Hall effect. Spatial profiles of spin accumulation and its evolution in an 
external magnetic field can be predicted by drift-diffusion modeling independently of any 
detailed treatment of the microscopic physics of the spin Hall effect. The practical ability 
to accurately model spin accumulation generated from the spin Hall effect enables reliable 
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engineering of electrically-generated spin accumulation profiles in devices. 
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Figure captions 

Figure 1. (a) Schematic of the 60-/xm wide channel with a transverse arm. The origin 
is defined to be at the center of the transverse arm in x and at the center of the main 
channel along y (red dot), and the dotted lines represent the location of scans discussed in 
the text, (b) 9k(B) taken away from the transverse arm for x 50/xm and y = — 26/xm 
(black circles), — 22/xm (black circles), and — 18/xm (magenta squares), (c) s z (y) for B = 
away from the transverse arm {x 50/xm). The inset shows the spin coherence time 
r(y) extracted by assuming a Lorentzian lineshape in the Hanle model, (d) Characteris- 
tic 9k(B) data from the center of the transverse arm at x — 0. The line is a guide to the eye. 

Figure 2. (a) Normalized 0k(B) from Fig. lc, corresponding to data near the channel 
edge, ~ L s from the edge, and ~ 2L S from the edge. Lines are calculations from the 
diffusion model described in the text, (b) r(E) and n{E) used in the model calculations. 
r is obtained from s z (B) scans near the channel edge for various while n{E) comes 
from measurements on a Hall bar. (c) Data is an image 9x{B,y) along the main channel 
width clearly showing the decreasing central peak width away from channel edge and sign 
change away from the central peak. The lower panel shows model calculations of s z (B ) y) 
as discussed in the text, (d) Data is an image 0K(B,y) at x = in the transverse arm 
showing the complicated spin accumulation profile. The lower panel shows diffusion model 
calculations of s z (B ) y) as discussed in the text. 

Figure3. Images of 6k and calculated s z (B = 0) in the transverse arm. The top panels 
are data and the bottom panels are model calculations as described in the text. The electric 
field in the channel has magnitude E = ll.lmV/ fim along x ) and is applied as (a) a square 
wave, (b) an offset square wave between and (c) an offset square wave between and 
—E. Data is normalized to the peak values for part (a). 
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Stern, ef al. 
Fig. 1 
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Stern, ef al. 
Fig. 2 
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Stern, etal. 
Fig. 3 



